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Abstract —Previous studies have confirmed the adverse im¬ 
pact of fading correlation on the mutual information (MI) of 
two-dimensional (2D) multiple-input multiple-output (MIMO) 
systems. More recently, the trend is to enhance the system 
performance by exploiting the channel’s degrees of freedom in the 
elevation, which necessitates the derivation and characterization 
of three-dimensional (3D) channels in the presence of spatial 
correlation. In this paper, an exact closed-form expression for 
the Spatial Correlation Function (SCF) is derived for 3D MIMO 
channels. This novel SCF is developed for a uniform linear array 
of antennas with nonisotroplc antenna patterns. The proposed 
method resorts to the spherical harmonic expansion (SHE) 
of plane waves and the trigonometric expansion of Legendre 
and associated Legendre polynomials. The resulting expression 
depends on the underlying arbitrary angular distributions and 
antenna patterns through the Lourier Series (LS) coefficients 
of power azimuth and elevation spectrums. The novelty of 
the proposed method lies in the SCL being valid for any 3D 
propagation environment. The developed SCL determines the 
covariance matrices at the transmitter and the receiver that form 
the Kronecker channel model. In order to quantify the effects of 
correlation on the system performance, the information-theoretic 
deterministic equivalents of the MI for the Kronecker model 
are utilized in both mono-user and multi-user cases. Numerical 
results validate the proposed analytical expressions and elucidate 
the dependence of the system performance on azimuth aud 
elevation angular spreads and antenna patterns. Some useful 
insights into the behaviour of MI as a function of downtilt 
angles are provided. The derived model will help evaluate the 
performance of correlated 3D MIMO channels in the future. 
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I. Introduction 

Multiple-input multiple-output (MIMO) systems have re¬ 
mained a subject of interest in wireless communications over 
the past decade due to the significant gains they offer in terms 
of spectral efficiency by exploiting the multipath richness of 
the channel. The increased spatial degrees of freedom not only 
provide diversity and interference cancellation gains but also 
help achieve a significant multiplexing gain by opening several 
parallel sub-channels. Pioneer work in this area by Telatar 
||T] and Foschini Q realized that capacity can potentially 
scale linearly with the minimum number of transmit (Tx) 
and receive (Rx) antennas for channel matrices with centered, 
independent and identically distributed (i.i.d) elements. These 
MIMO systems were designed to support antenna configu¬ 
rations capable of adaptation in the azimuth only. However 
recent measurement campaigns demonstrated that elevation has 
a significant impact on the system performance Q, Q. Ex¬ 
ploiting the channel’s degrees of freedom in the elevation can 
further enhance the system performance by benefiting from the 
richness of real channels. This has recently become a subject of 
interest among researchers and industrials. The reason can be 
attributed to its potential to open up possibilities for a variety 
of strategies like user specific elevation beamforming and cell¬ 
splitting. Encouraged by the initial implementations of this 
technology Q, the 3GPP is now working on defining future 
mobile communication standards in the frame of the study 
items on three-dimensional (3D) beamforming ||^. 

The discussion on the conspicuous advantages of 3D MIMO 
systems must be amalgamated with the observation that it is 
the orthogonality of the subchannels constituting the MIMO 
system that determines the extent of the multiplexing gain that 
can be realized. Large capacity gains can only be realized 
when the subchannels are potentially decorrelated. However 
in realistic propagation environments, the promised theoretical 
gains are not realized due to the significant spatial correlation 
present in the MIMO channel Q-GD- Therefore assuming the 
channel coefficients to be i.i.d is an oversimplification of the 
problems encountered in realistic propagation environments. 

The need to investigate the impact of spatial correlation 
on the performance of MIMO systems is acknowledged and 
well-known among researchers. However, most of the spatial 
correlation models proposed in literature are developed for 2D 
channel models and ignore the elevation ||^, ||T4|-||T8). These 
2D models assume an omnidirectional radiation of energy in 
the elevation. However, with the advent of smart antennas 
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GD, it is important to take into account the characteristics 
and patterns of the directional antennas that are being widely 
deployed now. These antenna patterns have often been de¬ 
scribed in literature using spherical harmonics | [20) , pT) , but 
this has not been done in the context of spatial correlation. A 
recent approach to incorporate these patterns in the correlation 
function was proposed in |22| , but the method could not 
admit a closed-form solution for this case and succumbed to 
numerical integration methods. 

The existing correlation models are derived for a particular 
distribution of the Angle of Departure (AoD) and Angle of 
Arrival (AoA) such as uniform, Gaussian, Von Mises, cosine 
or Laplacian 0, GD-Gl, @-|25). In fig, approximate 
closed-form expressions for the spatial correlation coefficients 
for clustered MIMO channel models were derived for Lapla¬ 
cian azimuth AoA distribution. The proposed method makes 
small angle spread approximation for uniform linear and circu¬ 
lar arrays and offers significant gains in terms of computational 
cost. Such assumptions on angular distributions can lead to 
useful closed-form expressions but do not accurately represent 
the characteristics of realistic propagation environments. In 
Hg, the authors derived exact closed-form expressions for the 
spatial correlation between Rx antenna elements for cosine, 
Gaussian and Von Mises azimuth AoA distributions. The use 
of Von Mises was shown to simplify the expressions and the 
impact of mutual coupling on the correlation was studied. 

The notion of spatial correlation in 3D propagation envi¬ 
ronments has been addressed in some research works. An 
important contribution in this area appears in | |2^ . The authors 
developed closed-form expressions for the spatial correlation 
and large system ergodic mutual information (MI) for a 3D 
cross-polarized channel model, assuming the angles to be 
distributed according to Von Mises distribution. The authors in 
showed that elevation plays a crucial role in determining 
the Spatial Correlation Function (SCF). The derivation is based 
on the spherical harmonic expansion (SHE) of plane waves and 
assumes the distribution of AoA to be 3D Von Mises-Fisher. 


In (251, closed-form expressions for spatial fading correlation 
functions of several omnidirectional antenna arrays in a 3D 
MIMO channel were derived and used for the evaluation of 
channel capacity. The derived results were expressed as a 
function of angular and array parameters and used to study 
the impact of azimuth and elevation angular spreads on the 
MI. However, this work assumes the angular distributions to 
be uniform. The analysis in f2S\ uses SHE of plane waves 
to compute the closed-form expressions for the correlation 
that can be applied to a variety of angular distributions. 
Although the tools presented are handy, the proposed closed- 
form solutions require certain assumptions to be made on 
the propagation environment. Even the simple assumption 
that the angles are uniformly distributed resulted in integrals 
involving Legendre polynomials that could not be expressed 
in a closed-form. Such assumptions neither aptly represent the 
characteristics of realistic propagation environments nor make 
the proposed method truly generic in nature. 

In order to quantify the effect of correlation, it is important 
to derive and simulate the correlated MIMO channels and 
characterize the information-theoretic MI for them. There are 


two widely used approaches to model these channels | [T6) . The 
hrst one is the parametric approach, in which the propagation 
paths are described using statistical parameters without being 
physically positioned. Channel realizations are generated by 
summing the contributions of multiple paths (plane waves), 
with specihc channel parameters like delay, amplitude, AoA 
and AoD. The second approach is nonparametric, wherein 
the SCE is used to determine the covariance matrices at the 
transmitter and receiver. These matrices are then employed to 
reproduce the spatial correlation across the MIMO channel. 
An example is the Kronecker channel model, which is useful 
for the evaluation of theoretical ML In this context, well- 
known results from Random Matrix Theory (RMT) have been 
employed to characterize the distribution of the MI of these 
channels in the asymptotic regime as the number of antennas 
at the base station (BS) and mobile station (MS) tend to 
grow large 0. Such theoretical results enable better 

understanding of the impact of the correlation on the MI. These 
so-called deterministic equivalents are reasonably tight even at 
moderate values of the number of antennas. 

The aims of this paper are fourfold. Eirst is to develop an ex¬ 
act closed-form expression for the SCE for 3D MIMO channels 
that can be used for any arbitrary choice of antenna patterns 
and distribution of azimuth and elevation AoD and AoA. The 
parametric 3D channel model used in the derivation is inspired 
from the models presented in standards like 3GPP SCM | [30) , 
WINNERh- | [3T| and ITU |(^. To get an analytically tractable 
closed-form solution, the SHE of plane waves and properties of 
Legendre and associated Legendre polynomials are exploited. 
The hnal expressions for the SCE developed for a uniform 
linear array of antennas are presented in Theorem 1 [@, 

(28 i]. The closed-form expressions depend on the underlying 
arbitrary angular distributions and antenna patterns through the 
Eourier Series (ES) coefficients of Power Azimuth Spectrum 
(PAS) and Power Elevation Spectrum (PES). To the best of 
authors’ knowledge, a SCE that works for the 3D channel 
model without making any assumptions on the underlying 
angular distributions and antenna patterns has not been de¬ 
veloped before. The second aim of this work is to validate 
the proposed SCE via simulations for angular distributions and 
antenna patterns specihed in the standards. The ES coefficients 
are computed and used to obtain the correlation coefficients 
that coincide with the Monte-Carlo simulated results. The 
third aim is to use the nonparametric Kronecker channel 
model for the evaluation of MI in the mono-user case. The 
developed SCE is used to determine the covariance matrices 
at the transmitter and receiver that form the Kronecker model. 
The pinhole phenomenon is discussed and illustrated as a 
restriction to the nonparametric Kronecker channel model as 
compared to the parametric channel model. The theoretical 
analysis for the mono-user case makes use of the deterministic 
equivalent of the MI presented in | [27) and studies the effect 
of angular parameters parametrized by azimuth and elevation 
angular spreads on the ML An interesting interplay between 
vertical antenna pattern and elevation spread is observed. 
The hnal goal of this work is to provide a Havor of the 
performance gains realizable through the meticulous selection 
of the transmit antenna downtilt angles in a multi-user scenario. 
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TABLE I. List of important symbols used in sections II-IV. 


Symbol 

Description 

e, I? 

Elevation AoD and AoA respectively. 

4>^ ‘F 

Azimuth AoD and AoA respectively. 

N 

Number of propagation paths. 


Complex amplitude of the path. 

^tilt 

Elevation angle of antenna Foresight. 


Tx antenna pattern. 


Rx antenna pattern. 


Horizontal antenna pattern. 

9t,V (^7 ) 

Vertical antenna pattern. 

kt 

Transmitted wave vector. 

kr 

Received wave vector. 

X 

Location vector of an antenna in 

V 

Direction of wave propagation. 

(t>SdB 

Horizontal 3 dB beamwidth. 

SsdB 

Vertical 3 dB beamwidth. 

Gp^max 

Antenna gain. 

Nbs 

Number of Tx antennas. 

Nms 

Number of Rx antennas. 

dt 

Distance between Tx antennas. 

dr 

Distance between Rx antennas. 

PAS 

Power azimuth spectrum. 

PES 

Power elevation spectrum. 

do 

Mean AoD in the elevation. 

a 

Angular spread in the elevation. 

Pt{s - s') 

Correlation between Tx antennas s, s'. 

Pr{u — u') 

Correlation between Rx antennas u, u'. 

jn 

Bessel function of order n. 

Pn 

Legendre polynomial function of order n. 

pm 

^71 

Associated Legendre polynomials. 

pm 


Pt, Pr 

‘^dt and '^dr respectively. 

Pnj p2ri,2k 

Legendre coefficients for even orders. 

j2m— 1 
^2n-1.2k-l 

Legendre coefficients for odd orders. 

^(f) ^<p 

ES coefficients of PAS. 

o-e,be 

PS coefficients of PES. 

No 

Number of terms summed over n 


in the proposed SCE 

P 

Mean AoD in the azimuth. 

K 

Angular spread in the azimuth. 


The MI analysis makes use of the deterministic equivalent 
of the signal-to-interference plus noise ratio (SINR) in | [33) 
with regularized zero forcing (RZF) precoding at the BS to 
mitigate inter-user interference. The researchers and industrials 
interested in using our correlation model need to provide only 
the FS coefficients of the PAS and PES they are using for the 
evaluation of their work. 

This paper is organized as follows. The 3D channel model, 
antenna configuration, PAS and PES are explained in Section 
II. In section III, we present an analytical derivation of the 
proposed closed-form expression for the generalized SCE. In 
Section IV, we provide simulation results that validate the 
developed SCE, adhering to most of the guidelines provided 
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in the standards. In section V, we present the nonparametric 
Kronecker channel model and recall well-known results on 
the deterministic equivalents of the MI of this model in the 
mono-user and multi-user systems. The performance of these 
systems is investigated as a function of channel and array 
parameters through numerical results. Einally in section VI, 
some concluding remarks are drawn. 

II. Channel Model and Power Spectrums 

Prior to proceeding into the derivation of the SCE for MIMO 
channels, it is vital to explain the characteristics of the 3D 
channel model and antenna configuration under investigation. 
MIMO systems of current LTE releases do not support antenna 
configurations capable of adaption in the elevation. However, 
encouraged by the potential of elevation beamforming to en¬ 
hance system performance, some standardized channel models 
have started to emerge that define the next generation 3D 
channels. We base the evaluation of our work on these channel 
models after making some realistic assumptions on the channel 
parameters. To improve the clarity of mathematical exposition, 
the important symbols used in sections II-IV are listed in Table 

m 

A. Standardized 3D Channel Model 

The MIMO channel model for which the SCE is derived 
is inspired from the standardized models like 3GPP SCM 
1^, ITU ||32) and WINNER These standards follow a 
system level, stochastic channel modeling approach wherein, 
the propagation paths are described using statistical parameters 
without being physically positioned. Channel realizations are 
generated by summing the contributions of multiple paths with 
specific parameters like delay, amplitude, AoA and AoD. 

Almost all system level based standards are 2D. However, 
owing to the growing interest in 3D beamforming, extensions 
of these standards to the 3D case have started to emerge 
recently in ||^, pT[ . Based on the aforementioned standards, 
the effective channel between the BS antenna s and the MS 
antenna u is given by |3T), |3g, 

N 

[H] 

su n \/9t{4^nt ^ni \lQr '^n)]u 

n—1 

^ ^n')'\s: (1) 


Azimuth plane 
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N 

[H]su = ^ an\/gt{(t>n,0n,0tiit)exp {ik{s - l)dtsm(j)nSm6n) y/ gr{(pn,'dn) exp {ik{u - l)dr sin (/?„ sini?„), (8) 

n—1 


where and are the azimuth and elevation AoD of 
the path respectively, (pn and are the azimuth and 
elevation AoA of the path respectively and OtM is the 
elevation angle of the antenna boresight. Note that Otm = 90° 
corresponds to zero electrical downtilt. is the complex 
amplitude of the path. The complex amplitudes are as¬ 
sumed to be i.i.d zero mean, ^ variance Gaussian RVs. Also 
yjgt{4>n,0n-,0tm) and are the global patterns 

of Tx and Rx antennas respectively where gt{(l>nj(^ntdtiit) ~ 
gt,H{4^n)gt,V and ~ gr,H i^n) gr,V 

Note that gt,H{(l>), gr,Hi‘f) '^^e horizontal antenna patterns 
and gt,v{d,dtiit),gr,vi^) are the vertical antenna patterns. 
Moreover, vectors at {(f>, 9) and ar(<p, '&) are the array responses 
of the Tx and Rx antennas respectively whose entries are given 

by, 

[at((^,6>)]s = exp(ikt.Xs), (2) 

[ar(v5,i?)]« = exp(ikr.x„), (3) 

where . is the scalar product, x^ and x^ are the location vectors 
of the Tx antenna and the Rx antenna respectively, kt 
and kr are the Tx and Rx wave vectors respectively, where 
k = ^v, with A being the carrier wavelength and v being the 
direction of wave propagation. Fig.[T]illustrates the 3D channel 
model being considered. It is evident that 9, Quit, "d G (0,7r) 
and 4),p & ’’■)■ 

B. Antenna Configuration 

If {er,ee,e^) is the spherical coordinate system, then verti¬ 
cal polarization refers to the polarization along eg and horizon¬ 
tal polarization refers to polarization along e^. Each antenna 
port comprises of vertically stacked antenna elements that 
determine the effective antenna port pattern. For the purpose of 
this work, vertically polarized antenna elements are considered. 
The antenna ports are placed at fixed positions along Cj,, with 
the elements in each port aligned along as shown in Fig. 
1^ There are Nbs and Nms antenna ports at the BS and MS 
respectively. The same Tx signal is fed to all elements in a port 
with corresponding weights to achieve the desired directivity. 
The MS sees each antenna port as a single antenna because all 
elements carry the same signal. Therefore, we are interested in 
the channel between the Tx antenna port and the Rx antenna 
port. 

In theory, the global antenna pattern of a port depends on 
the patterns of the elements within it and their corresponding 
weights. However to enable an abstraction of the role played 
by these elements to perform elevation beamforming, standards 
like 3GPP and ITU approximate the pattern of each antenna 
port by a narrow beam in the elevation. The combined pattern 
considered in dB is as follows [321 , [35) , 

Ap(fi, 9,9tiit) = Gp^max - min{-(A//((/)) -f Av{9, duu)), 20}, 

(4) 


where. 




Av{9,9tiit) 


—min 



20 


dB, 


—min 


12 


\ 9^dB 


2 

,20 


dB. 


(5) 


Therefore the horizontal and vertical antenna patterns at the 
transmitter can be approximated as, 

= -12 dB, (6) 

V 03dB J 

gt,v{9,9tat) = -l2(^-^r^\ dB. (7) 

Gp^raax = 17 dBl 'p 3 dB is the horizontal 3 dB beamwidth 
and 93dB is the vertical 3 dB beamwidth. The individual 
antenna radiation pattern at the MS, gripp'd), is taken to be 
0 dB in the standards. It is important to note here that we 
provided the expressions for the antenna patterns proposed 
by the standards only for the sake of completeness of the 
standardized 3D channel model. Also, they will be used to 
validate our proposed SCF later in section IV. However, the 
development of the SCF provided in this work is independent 
of the form of the underlying antenna patterns as long as 
gti^fint 9m 9tilt^ ~ gt,H{fin^gt^vi,9nj9^ill^ and gr{}Pn^’9ri) ~ 
gr,H {^n)gr,V (t^n)- 

Given the antenna configuration shown in Fig. v^.x^ = 
sin^„sin0„ and \r-Xu = sin(/j„ sini?„. The effective radio 
channel given in ([T]| can hence be written as (|^, where s = 
1, ■ • ■ Nbs and u = 1,... Nms- 



Fig. 2. Antenna configuration. 
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C. Power Azimuth and Elevation Spectrum 

Power azimuth spectrum (PAS) and power elevation spec¬ 
trum (PES) are important statistical properties of wireless 
channels and are shown to play an important role in determin¬ 
ing the spatial correlation present in the MIMO channel. They 
provide a measure of the power distribution upon the azimuth 
AoD and AoA and the elevation AoD and AoA respectively. 

Observing that the integral of the product of angular power 
density function of the azimuth AoD/AoA and the horizontal 
antenna pattern at the BS/MS yields the expected power 
transmitted/received by the directional antenna in the azimuth 
| [M1 , we define PAS at the transmitter as follows, 

PASi ((/))= ((/>), (9) 


III. Wavefield Decomposition and Spatial 
Correlation Function 

In this section, we rigorously derive a generic analytical 
expression for the SCF considering realistic antenna patterns 
and arbitrary distributions of AoDs and AoAs. 

Note that since the parameters describing the multipaths are 

1.1. d, the double sum in E[Hs„Hff^/] can be simplified to N 
multiplied by E[|ap] and the terms defined in (17i and (18i 
below. Therefore, it can be seen directly from (|^ that Iot 

1.1. d zero mean a’s with E[|ap] = the spatial correlation 
between the channels constituted by any pair of Tx and Rx 
antenna ports can be expressed as a product of the correlation 
between Tx antenna ports and the correlation between Rx 
antenna ports, i.e.. 


where the angular power density function p^{(j))=f^{(j>), the 
probability density function of azimuth angle. Therefore, 


P4>{(j))d(j) = 1. 


( 10 ) 


Similarly PES at the Tx side is defined as. 


PESt{e,euit) = gt,v{9,euit)pe{9). ( 11 ) 


The elevation angular power density function, Pe{ 9 )=^d^ 
which implies p7[. 




Pe{9) su\{9)d9 = 1. 


( 12 ) 


The same dehnitions and conditions can be extended to PAS^ 
and PESr. Note that the limits taken in ( |T2l i are (0, 27r) instead 
of (0, tt), which is the range over which 9 is dehned. This 
extension in limits, which would later assist in expressing 
SCF in terms of the FS coefficients of PES, entails that we 
dehne fe{9) to be zero from tt to 27r. This is generally 
true because the elevation angular density spectrums used in 
standards decay exponentially with 9. 

Example: The PES can be well fitted by the Laplace distri¬ 
bution p 8 |. The elevation angles are thus generated using. 


fe{9) oc exp 


V2\9-9o\ ^ 


sin 0 , 


(13) 


where a is the spread in the elevation direction and 9o is 
the mean AoD in the elevation. The density function decays 
exponentially and is zero for 9 ^ [ 0 , 7 r]. Hence to determine 
the constant of proportionality A, we can use condition ( [T^ 
and the observation that fg{9) = 0 for 0 ^ [ 0 , 7 r], 


SCF = E[[H«,h5„,] = pt{s - s')pr{u - u'), (16) 


where. 


pt{s - s') = E 
Pr{u — u') = E 


gt{(l),9,9tut)exp - s') sine/)sin 

(17) 

/ 27r 

gr{<p, d) exp li—dr{u — u') sin (p sin i/ 


(18) 


We derive the closed-form expression for the correlation 
between Tx antennas given in (171. This can be extended to 
the correlation between Rx antennas and the product of the 
two would yield the SCF between the channels. 


A. Spherical Hannonic Expansion of Plane Waves 

In a 3D propagation environment, the array responses of Tx 
and Rx antennas can be expanded using spherical decompo¬ 
sition for plane waves. Using the Jacobi-Anger expansion, a 
plane electromagnetic wave can be expressed as a superposi¬ 
tion of spherical waves p9] , 

OO 

g.fcx.v ^ ^ + l)j„(fc||x||)P„ (x.v) ,X e (19) 

n=0 


where k=^ is the wave number, v is a unit vector in the 
direction of wave propagation, x is the location vector of the 
antenna in jn is the spherical Bessel function of order n 
and Pn is the Legendre polynomial function of order n. 

We also state here the Legendre addition theorem of spher¬ 
ical harmonics 1 39 1 , | [4^ which will be employed later in the 
derivation of the SCF. When 7 is dehned as. 


Aexp An9d9 = 1, 


A = 


2 + a^ 


2\/2a sin9o + 2a’^e cosh 


(14) 

(15) 


cos 7 = cos 01 cos02 + sin 01 sin 02 cos((/)i — (/>2), (20) 

where ( 0 i, fi) and ( 02 , (/) 2 ) are the spherical coordinates of the 
vectors v and x respectively, then the Legendre polynomial of 
argument cos( 7 ) is given by ( [2T] i, where P™ are the associated 
Legendre polynomials. 
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n / _ \I 

P„(cos7) = P„(cos6»i)P„(cos6i2) + 2 V ^ ”^:^’ P™(cos6>i)P^(cos6*2) cos[w((/)i - ^ 2 )] 

' (n + m)\ 

m—l ^ ' 

( 27r \ cxD / 2 ^ \ 

exp f dt(s — 5 ^) sinc^sin^ j = ^ z^(2n + l)jn f 1^ ~ j i^n(sin(^sin^) 

Pt(s - s') =E\gti(l),O,0tiit)'^i""{2n + l)j„ ( -^dt\s - s'| ) (p„(cos 6 ')P„( 0 ) + 2 ^ ^ ^ ' 

L ^0 VA ;V + 

X P™(cos 6 »)P™( 0 ) cos I)) 


( 21 ) 

( 22 ) 

(23) 


Pt{s - s') = E [5t ((?!), 6 », 6 'ti;t)] jo (/3t|s - s'|) + ^(-l)”(4n + l)j 2 n {/3t\s - s'|) P2ni0)^P2nicOs9)gt,v{9,0tilt)]E[gt,Hi(l))] 

n—1 

oo / n \ 

+ ^ 4(-l)"j2n Wt\s - s'l) ^ (-l)™P|r(O)E[P|™(cos0)fft.v(0,0t,ii)]E[cos(2mj>)fft.^(</>)] (24) 

n—1 \ni—l / 

oo n 

+ ^ 4 z(- 1 )"j 2 „-i (/3t|s - s'l) ( ^ (-l)-P 2 “( 0 )E[P 2 “(cos%t,y(^ 9uit)nsm{{2m - ^gtM 


B. Spatial Correlation Function Using SHE of Plane Waves 

These results are now employed to derive a closed-form 
expression for the correlation between Tx antenna ports. The 
SHE result for plane waves in ( [T9] ) yields alternate expressions 
for the array responses of Tx and Rx antenna ports. It can be 
seen from Fig. 0 that for Tx antenna ports placed along 
direction with |]x|| = d(|s — s'|, x.v is given by sintjsinfl. 
Therefore the array response of the s*^ Tx antenna port can 
be expressed alternatively as ( |22| ). 

Also from Fig. the spherical coordinates {9i,4>i) of 
the wave vector k\ are (0, f), and the spherical coordinates 
{ 92 , 412 ) of X i.e., the vector along dt{s — s') are (f,f). 
Combining the addition theorem in (2T|i with (22), such that 
cos 7 =sin 4> sin 9 and using the resu ting expression in •ED 
would expand pt{s — s') to yield (23 1 . 

Proposition 1: For a uniform linear array of antenna ports 
with arbitrary antenna patterns and for arbitrary angular dis¬ 
tributions such that the 4>,(p £ [—tt, tt] and 9,-d G [0,7r], 
the correlation between any pair of Tx antennas ports can 
be expanded in a systematic way to yield (24i, where 

Pr(x)=^(n+i){gg;P™(x) and A = ^df 

The proof of Proposition 1 follows from the following 
properties of Legendre and associated Legendre polynomials, 

1) P„(0) = 0 if n is odd, 

2) P4f{0) = 0 if {n + m) is odd, 

3) Po{x) = 1. 

The following trigonometric relations were also used, 

1) cos(m'((j — |)) =(—1)™ cos{2m4>) for m' = 2m, 

2) cos(m'((j—^)) =(—1)"*“^ sin((2m—!)())) for m' = 2 to— 1, 
where m = 1,... ,n. 

Defining Pf^(x)=^{n + \) -Pr(x) and using the 

properties just described, (|23]) can be expanded in a systematic 


way. Finally, using the decomposition of gt(4>, 9,9uit) 


gt,Hi4>)9t,v{9,Stilt), and taking the deterministic terms out 
of the expectation operators yields ( |24| . 

The same approach would yield a similar expression for 
Pr{u — u') with gt{4>,9,9tiit) replaced by gr{ip,'9) and the 
AoDs replaced by AoAs. The expansion looks alarming but 
will be shown to yield an interesting closed-form expression. 


C. Closed-form Expression for SCF in terms of Fourier Series 
coefficients of PAS and PES 

The expansion in ( |24| ) exhibits several difficulties in deriving 
a closed-form expression for the SCF. The random variables, 
AoD and AoA, with respect to which the expectations need to 
be computed appear as the arguments of Legendre polynomi¬ 
als. Several tables of Legendre and associated Legendre poly¬ 
nomials exist that express the first few Legendre polynomials 
as functions of its arguments. However, we need a general 
representation that can be used for any order to facilitate the 
development of the expectation terms. For this purpose, we 
use the trigonometric expansion of Legendre polynomials pre¬ 
sented in ED- The following Lemma expresses the Legendre 
and re-normalized associated Legendre polynomials with even 
and odd orders as a linear combination of sines and cosines. 

Lemma 1 (From ED): For non-negative integers n and m, 

n 

P2n (cos X) =p^-\-2'^ Pn-kPn-\-k COs(2/i)x), 

k^l 

n 

Pin {cos x) = C 2 ™ 2 fc cos{2kx), (25) 

k^O 

n 

Pi:f_f\cOSx) =Y4nY2k-lM{2k-l)x), 


where the coefficients Pm c: 


2n 2k 


^ 2 n- 72 fc-i are generated 


using recursion relations in [|411, equations 2.8, 3.1-3.5]. 
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Pt{s - s') = 'K‘^a^{0)he{l)jo iPt\s - s'|) + ^(-l)”(4n + l)j 2 n {Pt\s - s'|) P2n{0)a^{0)Tr‘^ ^ Pn-kPn+k^[be{‘2,k + 1) 


k——n 


- be{2k - 1)] + ^ 4(-l)"j2n {^t\s - s'l) (-l)’"P|r(0)a^(2m)7r2 X] cl^, 2 k-M‘ik + 1) - 6,(2fc - 1)] (27) 

n—1 \m—1 k—0 / 

oo n n ^ 

+ x: 4 z(-i)"j2„-i ms - s'l) (E E 2 [««( 2 fc - 2 ) - ae{ 2 k)]] 




Pr{u - u') = 7r^a^(0)6,j(l)jo {Pr\u - u'\) + E(“l)"(4« + 1)72™ (/3r|w - w'|)P2n(0)a^(0)7r^ E Pn-kPn+k-^[bi)i2k + 1) 

n—1 k——n 

oo / ^ 1 \ 

- 6^(2fc - 1)] + E 4(-1)”J2„ {/3r\u - u'\) E (-l)'"i'lr(0)a^(2m)^2 E < 2^3 + 1) - 6,,(2fc - 1)] (28) 

n—1 \m—1 k—0 / 

00 n n .. 

+ E 4*(-l)”j2n-i m\u - ^'D ( E (-l)'"i^ 2 “( 0 )^( 2 m - 1)^2 ^ d“^_^-[a«(2fc - 2) - a^(2fc)]) 


k=l 


This Lemma is the most important ingredient in the deriva¬ 
tion of the SCF that gives it its generalized form that has not 
been derived before. We now state Theorem 1 that describes 
how the correlation between any pair of channels constituted 
by distinct pairs of Tx and Rx antenna ports can be computed. 

Theorem 1: For a uniform linear array of antenna ports with 
arbitrary antenna patterns and for arbitrary angular distribu¬ 
tions such that the (j),<p G [— 7 r, 7 r] and 6 , 1 } G [ 0 , 7 r], the SCF 
can be computed as. 


SCF = Pt{s — s')pr{u — u'), 


(26) 


where pt{s — s') and pr{u — u') are given by ( |27| ) and (28 1 
respectively, given that a^{k),b^{k),ag{k) and b^k) are the 
FS coefficients of PAS and PES respectively defined as. 


1 r 

aAm) = — / PASt(())) cos(TO()))d(^, (29) 

7^ l-TT 

1 r 

bAm) = — / PASt((/)) sin(m0)(i(/), (30) 

TI- J-K 

1 m 

ae{k) = - mSt{OAtm)cos{ke)de, (31) 

Ti- Jo 

1 m 

bg{k) = - PESt{e,9tM)sm{k9)de. (32) 

Ti- Jo 

The proof of Theorem 1 is postponed to Appendix A. The¬ 
orem 1 describes a novel method for obtaining the spatial 
correlation coefficients for 3D MIMO channels for arbitrary 
choices of antenna patterns and angular distributions which 
is often a difficult task. The proposed method is unique and 
is in contrast to most of the previous works that assume an 
underlying angular distribution and form of antenna patterns. 
This derivation can be generalized to other antenna topologies 
as well. 

Remark 1: The proposed SCF in Theorem 1 involves an 
infinite summation over n. However this infinite summation 


can be truncated to a small finite number, Nq, of terms such 
that the truncation error has a bound that decreases exponen¬ 
tially with Nq. This has been proved through an extensive 
analysis in | [42) and | |43) . The authors proved the bound for 
the 3D multipath field in (19 . We extend the analysis to the 
correlation expressions in (27i and ([28]l. 


pt{s - s') = E + ^)jn{k\\x\\)E[gt{(j),9,9tut)Pn (x.v)], 

n—0 

(33) 

^^0 = E *”( 2 ^+l)jn(^l|x||)E[gt(()', 6 >, 6 >tj/t)P„ (x.v)] 

n>No 


< E ( 2 ti + l)l 7 n(fc||x||)| |E[pt(()), 6 », 6 »tiit)P„ (sin())sin 6 »)] I 

n>No 

< E (2«+l)l7n(*l|x||)|Gp,,„a. (34) 

n>NQ 


since sup| 3 ,|<i |F’„(a:)| < 1 |44 , leading to 

\^[9t{(t>,9,9tiit)Pn{sin^sm9)]\ < Gp^rnax- @ is the 

same as equation (10b) in | |43| upto a scaling factor. The 
analysis in | [43) then uses the bound on spherical Bessel 
function and the Stirling bound on the Gamma function to 
show that for a finite ||x||, the multipath field and hence 
Pt{s — s') and pr{u — u') can be truncated to jnj < Nq, such 
that the truncation error is bounded by 


— G p^rnaxP exp(-(5), (35) 

where S = Nq — |"efc||x||/ 2 ], (5 > 0 and 77 « .678481234. For 
our analysis, to calculate the correlation between adjacent Tx 
or Rx antennas for an antenna spacing of ||x|| = .5A, given 
that Gp^max = 17dBi, Ao=14 would suffice to bound the error 
by approximately 0.5 %. 
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Fig. 3. Correlation between Tx antenna ports for distributions and patterns 
from standards. 


Fig. 5. Comparison of the coiTelation between Tx antenna ports for 2D and 
3D channels. 



Fig. 4. Correlation between Rx antenna ports for distributions and patterns 
from standards. 


IV. Validation of the SCF for Standardized 3D 
Channel Model 

In the last section, we presented a generalized methodology 
for obtaining a closed-form expression for the spatial correla¬ 
tion function for the 3D channel, that works for any arbitrary 
choice of antenna patterns and distributions of azimuth and 
elevation AoD and AoA. Since this paper largely focuses on 
the guidelines provided in the mobile communication standards 
used globally, it is important that our model is validated using 
the angular distributions and antenna patterns specified in the 
standards. 

In the standards, elevation AoD and AoA are drawn from 
Laplacian elevation density spectrum given in ( [T3] l, where 
9q is the mean AoD/AoA and a is the angular spread in 
the elevation. The characteristics of azimuth angles are well 
captured by Wrapped Gaussian (WG) density spectrum p2| , 
0, H). However in the recent years, the Von Mises 
(VM) distribution has received great attention in modeling 
nonisotropic propagation due to its close association with the 
WG spectrum p3) , | |45) . This distribution given by, 

exp(KCOs(a;-/r)) 

-■ (36) 


is related to the WG distribution through a straightforward 
relationship obtained using their first circulant moments p6t , 

WG(/i,cr^) = VM(/i,K), cr^ = 2[log/o(K) - log/i(K)], 

(37) 

where /„(k) is the modified Bessel function of order n, /i is 
the mean AoD/AoA and ^ is a measure of dispersion. 

In practice, directional antenna patterns are commonly used 
so the incorporation of these patterns in this work is note¬ 
worthy. As per the standards, the global antenna port pattern 
is given by © for each Tx port. However, in order to 
obtain a closed-form expression for the FS coefficients of 
PAS using VM distribution, the antenna ports are assumed 
to be omnidirectional in the azimuth, i.e. gt,H{(t>) = 0 dB. 
From a propagation viewpoint, it is the downtilt angle that 
determines the vertical pattern which would play a crucial 
role in determining the antenna directivity and its potential to 
change our perception on the physical distribution of scatterers. 
Therefore assuming non-isotropy in the azimuth will not affect 
the results to a considerable extent. The vertical antenna 
pattern gt,vi9,dtiit) is given by The pattern at the MS, 
gr{^,'9) = 0 dB in the standard^ecause mobile terminals 
should not, in most cases, favor any direction. Using these 
antenna patterns and angular densities, the FS coefficients of 
PAS and PES are computed and provided in Appendix B. 

The validation of the theoretical results in and ( [28] l 
is done by comparison with Monte-Carlo simulation results. 
The Monte-Carlo simulations are performed over two thousand 
channel realizations using (| 8 1 to compute the correlation values 
in and ([T8]). The angles of the propagation paths are 
generated using the densities described earlier and the antenna 
patterns from the standards are used for the purpose of 
validation. The Monte-Carlo simulated correlation between Tx 
antennas and Rx antennas is then compared with the theoretical 
correlation computed using ( [27| ) and ( [28| respectively. The 
theoretical results should provide a very accurate fit to the 
Monte-Carlo simulated correlation. For the simulations, we set 
No = 15, Otut = 95°, 03dB = 15°, hdB = 70°, a, = 7°, 
cr„ = 10°, 9o = 90°, Ks,Ku = 5 and fi = The results are 
shown in Fig. and Fig.|^ As expected, the correlation is seen 
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Fig. 6. Correlation between Tx antenna ports for uniform azimuth angular 
distribution and patterns from standards. 


to decrease as the distance between the pair of antenna ports 
increases and more importantly, it is apparent from the graphs 
that the theoretical results provide a perfect fit to the Monte- 
Carlo simulated correlation for as few as fifteen summations 
over n. 

A comparison is made against the 2D standardized model 
given by, 


[H]su = an^J gt,H{(l>n) exp {ik{s - l)dtsm(j)n) 

n —1 

\Jgr,H{.^n) exp {ik{u - 1 )dr sm ip n), (38) 

Note that this is obtained by using 0, d = 7 r /2 and 
gty{ 9 , 9 tiit), gr,v{d) = 0 dB in ([^. The correlation between 
the Tx antennas is then given by. 


pt{s - s') = E 


5 t,ff(<^)exp ( i^dt(s - s')sin(^ 


■ (39) 


The same theoretical analysis done for the 3D channel can 
be repeated here and would lead to the squaring of P 2 Tt( 0 ), 
P|™( 0 ) and and the removal of all expectation 

terms involving the elevation angles in ( [24| . The theoretical 
result for ( (39| ) is then obtained straightforwardly. The Monte- 
Carlo simulated correlation in ( [39l l and the theoretical result 
are plotted in Fig. The 2D model clearly overestimates the 
correlation, which is a consequence of ignoring the directivity 
of antennas and the propagation of multipaths in the elevation. 
Assuming the radiation of energy from all the antennas to 
be in the same fixed direction in the elevation will cause the 
antennas to appear more correlated. 

The developed results are further validated by using uniform 
angular distribution for azimuth angles. Both the horizontal 
and vertical antenna patterns, i.e. gt,H{(j>) and gt,v{d,9tiit) 
in and 0 respectively, are considered this time. The 
correlation obtained is real which results from the symmetry of 
the uniform distribution and horizontal antenna pattern causing 
& 0 (fc)’s to be equal to zero. Fig. illustrates the excellent 
agreement between our derived and Monte-Carlo simulated 
results and establishes the credibility of the proposed function. 


TABLE II. List of symbols in section V. 


Symbol 

Description 

Rbs 

Correlation matrix at BS. 

Rms 

Correlation matrix at MS. 

X 

Nms X ^BS matrix of i.i.d CAf{0, 1) entries. 

y, X 

Rx and Tx signals respectively. 

n 

Received noise (AWGN). 


Variance of n. 

K 

Number of users. 

G 

RZF precoding matrix. 

/3 

Scaling parameter of the precoder. 

C 

Regularization parameter of the precoder. 

gfe 

Precoding vector for the k*^ user. 

hfc 

Channel vector for the k^^ user. 


Data symbol for the user. 

Ik 

SINR for the user. 

Ik 

Deterministic equivalent of SINR. 

dLoS,k 

Elevation line of sight angle at the BS. 


V. Equivalent Channel Model and Mutual 
Information Analysis 

In this section, we focus on the MI analysis of the 3D chan¬ 
nel model in order to quantify the effect of fading correlation 
and measure the performance gains realizable through potential 
elevation beamforming at the transmitter by careful selection 
of the downtilt angles. The important symbols used in this 
section are listed in Table [III 

A. Kronecker Channel Model 

The idea that multiple antennas at transmitter and receiver 
can bring about remarkable improvements in the MI made 
MIMO methods exceedingly popular. However this improve¬ 
ment depends on the multipath richness since a large capacity 
gain can be realized in the presence of potentially decorrelated 
channel coefficients. It is therefore, important for channel 
models to take this correlation into account to allow for a 
more accurate performance analysis. There are two popular 
approaches to model these correlated MIMO channels. The 
first one is the parametric approach, which was discussed in 
detail in Section II-A. This approach takes into account the 
spatial characteristics of wireless channels quite meticulously. 
The second approach is nonparametric, wherein the spatial 
correlation in the MIMO channel is reproduced using theo¬ 
retical Tx and Rx spatial correlation matrices. The latter is 
more suitable for the information-theoretic analysis of MI. 
One of the widely used nonparametric channel models is the 
Kronecker model that relies on the two matrices describing the 
correlation characteristics at both ends of the communication 
link 0 , 0 , fig. This model is defined as, 

h = rLxrL. (40) 

where X is a NmsxNqs matrix whose entries are indepen¬ 
dently and identically distributed according to a complex cir¬ 
cularly symmetric Gaussian distribution, i.e. CV(0,1), Rms is 
the correlation matrix at the MS with [Rms]u,u' = Pr{u — u'), 
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^BS 


F((t^) = —— logdet I + —K{a^)RBs I + 2:7 — logdet I + —K(cr^)RMS )- 


Nbs 


(44) 


Rbs is the correlation matrix at the BS with [Rbs]s,s' = 
Pt{s — s') and pt{s — s'), pr{u — u') are obtained using the 
derived expressions in ( [Z7] l and ( [28) l. For antennas arranged in 
a linear array, Rbs and Rms are Toeplitz. 

We consider the downlink of a single cell, where the BS is 
equipped with Nbs antenna ports and the MS is equipped with 
Nms antenna ports. The channel is linear and time-invariant. 
A time-division duplex (TDD) protocol is considered where 
BS acquires instantaneous CSI in the uplink and uses it for 
downlink transmission by exploiting the channel reciprocity. 
The channel H is known only to the receiver but not to the 
transmitter. Therefore power is distributed equally over all 
Tx antennas instead of employing the water-filling scheme. 
Moreover H is fixed during the communication interval, so 
we do not need to time average the MI. The received complex 
baseband signal y S ]y[§ gjyen by. 


y = Hx-t-n, (41) 

where x G is the Tx signal from the BS, H is the 

iVMSxiVss channel matrix generated using (40i and 
is the additive white Gaussian noise (AWGN) with variance cr^. 
The MI of the Nms'>^Nbs MIMO system with equal power- 
allocation is then given by. 


/(cr^) = log det Inms + 


1 


Nbs 


HH 


H 


(42) 


where cr^ is the noise variance and the average total Tx power 
is assumed to be 1. 

The subsequent analysis in the next two subsections will 
make use of the deterministic equivalents of the MI available 
in literature to study the behavior of the correlated 3D MIMO 
channels through numerical results. In the numerical results 
that follow hereafter, we set Nbs = 20, Nms = 20, 
Wo = 15, 9tm = 95°, O^dB = 15°, hdB = 70°, a, = 3°, 
ffu = 10°, 00 = 90°, Kg = 5, Ku = 5 and p = 0. For 
the mono-user systems, we work for the SNR level of 0 dB. 
Two thousand independent Monte-Carlo realizations of the 
parametric channel from ([^ are generated to compute the MI 
values in ( |42j l to allow for comparison with the theoretical 
results that follow. For these simulations, the azimuth angles 
are again drawn from VM distribution given in ( [36] l and the 
elevation angles are drawn from Laplacian distribution given 
in (13 I. The vertical antenna pattern given in Q is considered 
but again the antennas are assumed to be omnidirectional in 
the azimuth to allow for the calculation of FS coefficients of 
PAS in a closed-form. The theoretical correlation coefficients 


are calculated using (27 1 and (281. The developed SCF is used 
to determine the covariance matrices at the transmitter and 
receiver which are needed for the Kronecker model in (|40ll. 


B. Mutual Information Analysis of a Mono-User System 

It is imperative to study the behaviour of MI of MIMO 
channels in the presence of fading correlation to evaluate 



Fig. 7. Pinhole effect. 


different beamforming techniques. The MI for every realization 
of the channel can be viewed as a RV and it is interesting 
to study the statistics and distribution of this RV. Deriving 
closed-form expressions for the distribution of MI of the 
Kronecker model is a challenging task. However in the large 
(Nbs, Nms) regime, RMT provides some simple determin¬ 
istic approximations to this distribution. These deterministic 
equivalents are quite accurate even for a moderate number of 
antennas. The deterministic equivalent for the MI of Kronecker 






model was studied by |^, where it was shown that 
converges to a deterministic quantity defined as the fixed 
point of an integral equation. More recently, Hachem et al 
derived a deterministic equivalent for the MI of Kronecker 
channel model and rigorously proved that the MI converges to 
a standard Gaussian random variable in the asymptotic limit in 
The result is stated here and will be used in our analysis. 

Theorem 1 in p7| . Theorem 4.1. in | [28) suggests that for 

the Kronecker channel model given by H=R|,jgXR^g, the 
ergodic MI converges to a deterministic quantity, given that 
the assumptions Al and A2 in p7) on the eigenvalues of R^s 
and Rms hold. Under this setting, the authors showed. 


Nbs 


E[Iia^)]-Via^) 




-^ 0 , 


(43) 


where V (cr^) is given by (WS and where, (k((T^), K{a^)) is the 
unique positive solution orme system of equations given by. 


K{a^) = 


Nbs 


-tr R 


MS 


I-f —K{a^)RMS 


(45) 


iiia^) = 


Nbs 


tr I Rbs (I + ^K(cr^)RBS 


(46) 


This theorem shows how the MI can be approximated by a 
deterministic equivalent in the asymptotic regime. The result 
remains relevant even for a moderate number of antennas as 
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Normalized antenna spacing, dA 


Fig. 8. Verification of deterministic equivalent of MI of Kronecker model. pjg j j elevation angular spread on correlation. 


-t(- Monte-Carlo, Kj= 5 
-e- Theoretical, k^=5 
Monte-Carlo, Kj= 30 ~ 
-^Theoretical, Kj= 30 
Monte-Carlo, Kj= 45 
-^Theoretical, k=45 



0.5 1 1.5 

Normalized antenna spacing, d^/A, 


Fig. 9. Effect of azimuth angular spread on correlation. 



Fig. 10. Effect of azimuth angular spread on mutual information. 


will be confirmed through simulations. Before discussing the 
results, it is crucial to point out a restriction to Kronecker 
channel model regarding the loss of information about the 
number of paths. The rank structure of MIMO channel matrix 
not only depends on the correlation present in the channel 
but also exhibits a strong dependence on the structure of 


scattering in the propagation environment. It is possible to 
have a rank dehcient channel matrix even if the fading is 
decorrelated at both ends due to a small number of multipaths 
as compared to the number of antennas. This phenomenon, 
known as pinhole or keyhole effect, is generally observed in 
mild scattering conditions or when the communication link is 
very long. It has been studied in pT) , pS) . The pinhole effect 
is captured by the parametric channel models, like the one in 
that explicitly depend on the number of propagation paths. 
However the nonparametric model in (j4^ does not exhibit this 
phenomenon as there is no information about the number of 
multipaths involved. This has been illustrated in Fig. Monte- 
Carlo simulations are performed over two thousand channel 
realizations using and ( [40| to compute the MI values for 
the parametric channel model and non-parametric Kronecker 
channel model respectively for different number of propagation 
paths N. It can be seen that a reduction in the number 
of paths can severely affect the channel MI as conhrmed 
by the results for the parametric model. However, the MI 
obtained using the Monte-Carlo realizations of the Kronecker 
channel model is unaffected and only depends on the degree 
of spatial correlation present in the channel. Next we verify 
the deterministic equivalent for the MI of Kronecker channel 
model in ( [43] l. The theoretical result in ( |44| i is plotted in Fig. 

for TV = 20 and is seen to coincide quite perfectly with the 
MI obtained using the Monte-Carlo realizations of the non¬ 
parametric Kronecker model in (401, for a moderate number 
of antennas. The MI using the Monte-Carlo realizations of the 
parametric model in ([^ is also plotted and the pinhole effect 
is again apparent. 

We now use the deterministic equivalent of MI to study the 
impact of angular spreads on the MI. It has been observed 
in the past that higher azimuth spread can cause the adjacent 
antennas to appear close to uncorrelated. This has been con¬ 
hrmed in the results shown in Fig. and Fig. 10 An increase 
in K corresponds to a decrease in the spread which causes the 
correlation to increase. A decrease in the channel MI is then 
evident. Only the normalized theoretical MI in (44i for the 
nonparametric model has been plotted for = 40 to minimize 
the pinhole effect. The number of antennas at the BS and the 
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Fig. 12. Effect of elevation angular spread on mutual information. 



Fig. 13. Effect of elevation angular spread on MI for users away from the 
direction of antenna boresight. 


MS are increased at the same rate as per Telatar’s finding that 
the MI scales with minCTVss, Nms) 0- Since the correlation 
function is not insensitive to elevation, so in the context of 3D 
channels it is important to analyze the behaviour of MI with 
changes in angular parameters in the elevation. The correlation 
still decreases with an increase in the elevation angular spread 
as shown in Fig. 11 However the result for MI reveals an 
interesting interplay between the Tx power and the spatial 
correlation which is a consequence of the incorporation of the 
antenna pattern in our model. Note that the value of spatial 
correlation at 0 antenna spacing is nothing but the average 
Tx power of the MIMO system. It can be seen from Fig. 


11 that an increase in the elevation spread can undoubtedly 


cause the correlation to decrease. However, for users in good 
conditions with Oq « Ouit, i-e. for users within the direction of 
the antenna boresight, the incorporation of the antenna pattern 
into our channel model causes the Tx power to reduce with 
the increase in the spread as seen through the values of pt 
at dt/X = 0. This results in an overall decrease in the MI as 
shown in Fig. [T^ This effect has not been observed in previous 
works that deal with the elevation because they generally do 
not consider antenna patterns in their results for SCFs and 
MI. It is important to note here that the incorporation of the 
antenna pattern would not change the expected behaviour of 
the increase in MI with the increase in the spread for users in 
bad conditions, i.e. where the user is located far away from the 
elevation angle of the antenna boresight. In that case, a higher 
spread would actually benefit this user as it will manage to 
receive the signal from some propagation paths with significant 
energy and will benefit from the higher spread. This has been 
illustrated in Fig. 13 for the user receiving energy from the 
mean AoD 6 q = 130° while Quu is still equal to 95°. 

Finally we use the theoretical equivalent of MI in (|Hi to 
compare the MI of the 2D and 3D Kronecker channel models. 
The former is obtained by using the theoretical solution of ( 3^ 
to obtain the entries of [R_bs]s,s' = pt{s — s'). Similarly, the 
entries of [Rms]ii,«' are obtained by getting the theoretical so¬ 
lution of Pr{u — u') = E [gr.Hi^f) exp {i'^dr{u — u') sinyi)1 , 
after assuming d = 7r/2 and gr.vid) = 0 dB in The 
result is provided in Fig. and confirms that incorporating 



Fig. 14. Comparison of the MI of the 2D and 3D Kronecker channel models. 


the elevation affects the estimated values of the MI to a 
considerable degree. The lower MI for the 3D channel model 
is explained by the lower amount of transmitted power in the 
3D case as seen through the value of pt at dt/X — 0 in Fig. 
1^ This phenomenon is a result of incorporating the antenna 
pattern in the 3D model and is observed for users located in 
the direction of the antenna boresight as explained earlier in 
the discussion on the effect of angular spreads on the MI. 


C. Mutual Information Analysis of a Multi-User System 

A more robust system to channel correlation is the multiple 
user MIMO system, wherein instead of using multiple antennas 
for a single receiver, multiple users are served simultaneously. 
An obvious disadvantage of this system is the inter-user inter¬ 
ference which necessitates the use of an intelligent precoding 
scheme at the transmitter. Regularized zero-forcing (RZF) is 
one of the state-of-the-art schemes that we will employ to 
mitigate inter-user interference in our analysis of multi-user 
system. The RZF precoding matrix at the BS is given by, 

G = (HH^ + CNsslssy' H, 


(47) 
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a* AA fifi A 

K users 


Fig. 15. Multi-user scenario. 



Fig. 16. Effect of tilt angle on mutual information in a multi-user scenario. 


where is a strictly positive regularization parameter, /3 is a 
scaling parameter such that tr(G^G) — P, where P is the 
total available Tx power and K is the number of non co¬ 
operating users such that Nbs > K to avoid user scheduling. 
The channel vector for user G is given by, 

hfc = (48) 

where has i.i.d zero mean, unit variance complex Gaussian 
entries, RBS,k is the per user Tx correlation matrix and, 

Ptx X PLfe X AG X SF 

Qk = -s-, (49) 

cr^ 

where PL^ is the path loss experience by user k, SF is the 
shadow fading, Ptx is the transmitted power and AG is the 
antenna gain. PL^ is computed using the path loss model 
proposed for Urban Macro (UMa) scenario in 1^, SF =6 dB, 
antenna gain=17 dBi and cr^=1.13 x ]32|, | |4^ . 

The BS uses linear precoding. The precoding vector for 
the user is G and the data symbol is 

Sk CAf(0,l). Therefore the BS transmits the Nbs x 1 signal. 


K 

^ = J2^kSk = Gs. 


(50) 


The SINR for the k*^ user is then given by. 


hfgfcgfhfc 


Ik = 


hfGG"hfc-hfgfcg"h, + l 






(51) 


In this setting, a deterministic equivalent can be computed 
for the SINR, 7 ^ for every user k using tools from RMT, such 
that the convergence to the deterministic equivalent is almost 
sure as the number of antennas at BS and K tend to infinity 
p?j. T his deterministic equivalent was derived in great detail 
in 11331 (Theorem II). The theorem says that. 


7fc - 7fe 


-^ 0 

Nbs 


(52) 


where 7 ^ is defined in [|^, equation (19) using equations 
(20)-(25)]. Therefore the deterministic equivalent of the MI of 
user k in the multi-user case would be given by. 


/(a) = log(l + 7D- (53) 


This convergence result holds even for a moderate number 
of antennas and users and facilitates the analysis of the impact 
of downtilt angles on the system performance in a multi-user 
system. We analyze the multi-user scenario with iT = 40 and 
Nbs = 60. The scenario is illustrated in Fig. 15 The users 
are randomly positioned between radii of 100m and 250m 
from the BS. Qk is computed using ( |49l l for every user and 
the deterministic equivalent of SINR is computed based on 
Theorem II in p3j . This helps us get a theoretical equivalent 
for the Ml for every user. The elevation line of sight (LoS) 
angle with respect to the horizontal at the BS is computed for 


each user using OboS^u = tan 


-1 


Ati 


=, where {xk,yk) are 


the coordinates of the user k and A/i is the height difference 
between MS and BS. The mean elevation AoD used in the 
computation of FS coefficients of PES is also equal to dLoS,k- 
The users are located with dLoS,k G [95.37°, 103°]. We piot 
in Fig. 16 the normalized Ml (Total MI of the system/AT) of 
this system using the theoretical deterministic equivalent from 
and also using the Monte-Carlo realizations of SINR in 
and see that the deterministic equivalent is quite accurate 


(53 


(51 


even for a moderate number of users and BS antennas. More 
importantly, the result shows that the performance of the users 
at the cell edge is most sensitive to the value of downtilt angle. 
The performance of the system is maximized when the antenna 
Foresight angles at the BS are set equal to the OloS of the 
user at the cell edge, i.e. Btm ~ 96°. This result highlights 
the prospects of elevation beamforming in enhancing system 
performance of the future correlated MIMO systems. 


VI. Conclusion 

In this paper, we characterized the spatial correlation present 
in a 3D MIMO channel, assuming a uniform linear array of 
antennas. The conventional SCFs do not take into account 
the effect of elevation and antenna patterns, which renders 
them unsuitable for the evaluation of future correlated 3D 
MIMO channels that are currently being outlined in the 
next generation of standards. We derived the proposed SCF 
using SHE of plane waves and properties of Legendre and 
associated Legendre polynomials. The final expressions for the 
SCE presented in Theorem 1 illustrate how this generalized 
function depends on the underlying arbitrary antenna patterns 
and angular densities through the ES coefficients of PAS and 
PES. Numerical results show an excellent agreement between 























IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. XX, NO. X, XX 2014 


14 


the derived theoretical and Monte-Carlo simulated results for 
the spatial correlation. Furthermore, to quantify the effects 
of correlation on the system performance, we study the MI 
of the nonparametric Kronecker channel model. This model 
conhrms the existence of pinhole channels and allows us to 
use the available deterministic equivalents of the MI in the 
asymptotic limit for both mono-user and multi-user cases. By 
expressing the SCF in a closed-form as a function of channel 
and array parameters, the impact of azimuth and elevation 
angular spreads on the MI is investigated. The simulation 
results also provide useful insights into the impact of antenna 
patterns and antenna tilt angles on the achievable rates and 
conhrm the potential of elevation beamforming to enhance the 
system performance. 


Appendix A 
Proof of Theorem 1 

To derive an analytical expression for the SCF, we use the 
trigonometric expansion of the Legendre polynomial (25 1 |411. 
To this end, note that the expectation terms involving elevation 
angles in (24i can be written using the condition in (12i as, 

p27T 

E[P 2 nicos 9)gty{0,6uit)] = ^ ^pn-kPn+k J cos(2A:6») 


k——n 


X gt,v{9,9tiit)p{9) sin{9)d0 


= E 


Pn—kPn-\-k 


k——n 
1 /•27r 




sm{{2k + l)e)PESt{e,etm)d9 


sini{2k - l)0)PESti0,0tM)d9 


(54) 


E[P|™(cos(0))5t.y(0, 0tut)] = E [< 2 k / cos(2fc0) 




X gt,v{0, 9uit)p{0) sin(6»)d0 


= E< 


,2 m 
'2n,2k 


r27T 


sm{{2k + l)0)PESt{0,9tm)d0 


1-27! 


sin{{2k-l)0)PESt{0,0uH)d0 


(55) 


n 

E[P^::ri\cos{9))gt,v{0,0uit)] = E [^“. 2^-1 




X / sm{{2k-l)e)gty{0,6tiit)p{0)sm{6)d0 

Jo 

1 p27T 

= E 4n2k-i [2 / - 2)0)PES*(0, 9uit)d9 

1 _1 0 


cos{{2k)0)PESt{0,0tm)d0 


(56) 


As explained in section II-C, the limits in (|54)i-(|56]l have been 
set to (0, 27r) as fg{0) is 0 outside (0, tt). Doing this, we have 


expressed the expectations involving the elevation angles as a 
linear combination of the scaled FS coefficients of PES. Eor 
the expectations involving azimuth angles. 


E[cos(2m0)5t,ff((/')] = [ cos{2m(l))gt,H{4’)f{4’)d(j) 

J —TT 

cos(2TO^)PASt(^)d(/). 


(57) 


E[sin((2m - l)(j))gt,H{(t))] 
X f{(j))d(j) = f sin((2TO 


sin((2m 




l)(^)PASt((^)ci</.. (58) 


They express the expectations involving azimuth angles as a 
linear combination of scaled ES coefficients of PAS. Defining 
the PS coefficients of PAS and PES as. 


1 r 

aAm) = — / PASt{(j)) cos{m(j))d(j), (59) 

^ J-7V 

1 

bPm) = — / PASt(0) sin(TO(^)(i(^, (60) 

7^ J -77 

1 

ae{k) = - j PESt{0,0uit)cos{k0)d0, (61) 

TT Jo 

1 

beik) = - PESt{9,9uit)sm{k9)d9, (62) 

77 Jo 


one can immediately see the relationship between the terms 
involving expectations and these coefficients. Plugging ( [59|- 
(j6^ in and using the resulting expressions in ( |24| i 

would result in a closed-form expression for pt{s — s') yielding 
(27 1 . A similar development would yield (2^ by replacing 
AoD with AoA and gt replaced with This completes the 
proof of Theorem 1. 


Appendix B 

PouRiER Series Coefficients of PAS and PES 


A. FS Coefficients of PES at BS 

The PS coefficients are computed assuming Gp^max = 0 
dB. The value of Gp^max = 17 dBi will be incorporated in 
the simulations as a scaling factor (Antenna Gain). 


I \ ^ r ( 

ag,i{m) = — exp - 
77 Joo \ 


V2{9-0o)' 


X exp I —1.2 ( log(lO) ) cos{m0)d9, 


A 


3? 


/avT 
X ( erf 


exp —c 


2 Trajm + b'^ + 2bjm — mj 


2a| + b + jm 

2 y/d 


- erf 


4a 

2a(| — 0o) -\-h + jm 
2i/a 
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where a = b = log(io) ^ ^ ^ ^ 


1.2(f-^t,zt)2log(10) V2{^-9o) 


^idB 


^idB 


and erf is error function. 


,, A ( 72(00-0)' 

doA^) = — - 


X exp 1.2 

A 


3? 


aiT 


)_^ \ \ 

——— ) log(lO) I cos{m9)d9, 

bsdB ) j 

2 'Kajm + 5^ + 2 bjm — 


exp —c 


X erf 


2 a(| - 0 o) + b + jm 

2y/a 


4a 


- erf 


—2a| + b + jm 
2yAl 


be expressed as a series of Bessel functions as, 


PAS{(j)) = ^ 1^1 + ^ Ij{k) cos(j(«i' - Ai)) j , (65) 

1 ^ oo 

= — + , ' V Ij{K)[cos{j(j)) cos(j» + sin(j(/)) sin(j»] 

iTT TT1o[K) 

J = 1 


a^{m) = ImiK)cos{mfj.), 

7r/o(K) 

b<p{'m) = . . I^{K)sm{mn). 


( 66 ) 

(67) 

( 68 ) 


where a = 

^3dB 

1.2(f-gt.i,)^log(10) 


b = 


-2.4(f-et,it)log(10) 


V2(f —So) 


^idB 


— — and c = 

(7 


Also 5ep(m) and 6^,2 (w) have the same expressions as 
a 8 ,i{m) and ae^ 2 i'iTi) respectively with only 5ft replaced with 
Sr. The FS coefficients are. 


asim) = a 0 ^i{m) + ag^ 2 {m), (63) 

bg{m) = b 0 ^i{m) + bg^ 2 {m). (64) 


B. FS Coefficients of PES at MS 
Note that griff) = 1- 


asim) = — [ exp [ — 

Jo \ 




cos (m0) d9, 


r 272 , . , f -tt 

= M2 + mV0 *’»«.) - V“P (vS 

X (exp + (-ir exp (=^7^ 


b^im) = — / exp ( - 

Jo \ 


72|t^-0ol 


Aa'^ [272 . . „ , 

sin im9o) + m exp 


sin (to0) d9 

( 


tt{ 2 + m‘^a‘^) I a 
'72(1- 0 o) 


X exp 


-(-L'^exp 


V72ax 

-72(1- 0 o) 


C. FS Coefficients of PAS 

Since the antennas are considered omnidirectional in the 
azimuth, so gt^nif) and 3 r,H(‘p)=l- The VM distribution can 
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